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Abstract. In this paper, we use a traveling wave reduction or a so-called spa- 
tial approximation to comprehensively investigate periodic and solitary wave 
solutions of the modified Benjamin, Bona & Mahony equation (BBM) to in- 
clude both dissipative and dispersive effects of viscous boundary layers. Under 
certain circumstances that depend on the traveling wave velocity, classes of 
periodic and solitary wave like solutions are obtained in terms of Jacobi ellip- 
tic functions. An ad-hoc theory based on the dissipative term is presented, in 
which we have found a set of solutions in terms of an implicit function. Using 



dynamical systems theory we prove that the solutions of (1.6 1 experience a 
transcritical bifurcation for a certain velocity of the traveling wave. Finally, 
we present qualitative numerical results. 
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1. Introduction 

The first published observation about a solitary wave propagating along a uniform 
canal was made by John Scott Russel in 1834, who observed a solitary wave with 
an amplitude of about a foot and a width of about thirty feet on the Edinburgh to 
Glasgow canal. In his report to the British Association [T], Russel describes how he 
followed the solitary wave for more than one mile on horseback, observing that the 
observed solitary wave displayed a remarkable property of permanence preserving 
its original shape. He also noticed that taller waves will travel faster, and that 
solitary waves that move with different speeds, undergo a nonlinear interaction 
from which they emerge in their original shape. 
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In 1872, Joseph Boussinesq proposed a variety of possible models for de- 
scribing the propagation of water waves in shallow channels \2\ , including what is 
referred to as the Korteweg de Vries (KdV) equation 

In 1895, Korteweg and de Vries, under the assumption of small wave ampli- 
tude and large wavelength of inviscid and incompressible fluids, derived an equa- 
tion for the water waves, now known as the KdV equation, which also serves as 
a justifiable model for long waves in a wide class of nonlinear dispersive systems. 
KdV it has been also used to account adequately for observable phenomena such 
as the interaction of solitary waves and dissipationless undular shocks. For the 



water wave problem, (1.1 1 is nondimensionalized, since the physical parameters 
|u| = ?>rj/2H, X = \f^x* jH, and t = ^J&g/Ht* ^ where rj is the vertical displace- 
ment of the free surface, H is the depth of the undisturbed water, x* is dimensional 
distance and t* the dimensional time are all scaled into the definition of nondimen- 
sional space a;, time and water velocity u{x,t). When the physical parameters 
and scaling factors are appropriately absorbed into the definitions of u, x and t, 
the KdV equation is obtained in the tidy form 

(1.1) 



A further reduction could be made by removing the second term of (1.1) by tak- 



ing x' — X — t and t as independent variables, but nothing of significance is ac- 



complished by this. Eq. (1.1 1, or its equivalent without the second order term, 
is commonly taken as the staring point for mathematical studies of long-wave 
phenomena, although facts with considerable theoretical significance are already 
entailed in the derivation of KdV, [3]. Thus, the condensed form tends to disguise 
the meaning of the theory of the equation with regard to the original physical 
problem. 



Although (1.1 1 has some remarkable properties it manifests non-physical 



properties; the most noticeable being unbounded dispersion relation. It is helpful 



to recognize that the main difhculties presented by ( 1.1 1 arise from the dispersion 
term and arise in the linearized form 

Ut +Ux: + Uxxx = (1.2) 



First note that when the solution of ( 1.2 ) is expressible as a summation of Fourier 



components in the form F{k)e ^ the dispersion relation is 

uj = k-k^. (1.3) 

The phase velocity ui/k becomes negative for k'^ > 1, in contradiction of the 
original assumption of the forward traveling waves. More significantly, the group 
velocity ^ = 1 ^ S/c^ has no lower bound. 

To circumvent this feature, it has been shown by Benjamin, Bona & Mahony 



[1] that (1.1 1 has an alternative format, which is called the regularized long- wave 
or BBM equation 

Ut + Ux + uux ~ Uxxt ^ 0, (1-4) 
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in which the dispersion term u^xx is replaced by —Uxxt that results in a bounded 
dispersion relation. That was utilized to prove existence, uniqueness, and regularity 
results '3]. It is contended that (1.4 1 is in important respects the preferable model 
over which is unsuitably posed model for long waves. They also showed that 
(1.4) has the same formal justification and possesses similar properties as that of 
KdV and its solutions generally have better smoothness properties that those of 
(1.1). The linearized version of (1.4) has the dispersion relation 

^ (1.5) 



l + fc2^ 



according to which both the phase velocity uj/k and the group velocity ^ are 
bounded for all k. Moreover, both velocities approach zero for large k, which implies 
that fine scale features of the solution tend not to propagate. The preference of 
(1.4) over (1.1) became clear in [1], when the authors attempted to formulate an 
existence theory for ( 1.1 ), respective to nonperiodic initial condition u(x, 0) defined 
on (— oo, oo). In contrast, the existence and stability theory for (1.4) is essentially 
straight forward and perhaps the most persuasive evidence that (1.4 1 is better 
founded than (1.1) as a convenient model. 

A generalization to ( 1.4 ) to include a viscous term is provided by the equation 

Ut+Ux+ UUx - Uxxt = VUxx (1-6) 

where v is transformed kinematic viscosity coefficient of a liquid. Here, we will 
study the above equation and henceforth we shall refer to (1.6) as the modified 
BBM. 

The structure of the paper is as follows. In §2 a traveling wave solution of 
(1.6) will be derived. Depending on the traveling wave velocity c, three types of 
analytical solutions are found. When the viscosity is present a solution in implicit 
form is found using an ad-hoc theory for which the results will be shown numer- 
ically. In §3 we will closely look at the principles of linearized stability for the 
traveling wave solution using dynamical systems theory for the modified BBM. 
We will conclude the paper with the discussions and conclusions in §4. 



2. Traveling solutions to the modified BBM 

The class of solitons solutions is found by employing the form of the traveling wave 
solutions of (1.6) which takes the form of the Ansatz 

w(x,t) = 0(C) (2.1) 

where C, = x — ctis the traveling wave variable, and c is a non negative translational 
wave velocity. The case when c < 0, i.e., solitons traveling in the opposite direction, 
can be treated in a similar fashion by letting c — >■ — c. The substitution of (2.1 1 in 
(1.6) leads, after some simplification to 

(1 - c)0c + ^(0')c + - ^-^CC = (2.2) 
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Assuming that 0, 0^ — > as |C| — ^ ±oo, and by integrating once (2.2 1 we obtain 
the Lienard equation 

20CC = «'^C - Z?-^' - 7-^, (2.3) 

where a = — , /3 = -, and 7 = , which can be put in the first order form by 

making the substitution Lo{<p) — 4>i^ 



(2.4) 



Then, by making the substitution ri{(f)) — ^a;(0), the last equation may be written 
in the form 

W0 - ?7 = + 0.41, (2.5) 

where a = — and & = — We recognize this as the Abel's equation of the 
second kind written in the canonical form. 



2.1. No viscosity, ly = 

If the viscosity is not present, i.e., 1^ — 0, then a — 0, and hence eq. (2.2) becomes 

1 



2c 
1 + c 



c- 1 

: ( 

C 

: — Cd) 



if 



if 



c > 1 
< c < 1 



with analytical traveling waves solutions 

u{x,t) =3(c- l)Sech^ 
3c 



u{x, t) 



1 + c 



-Sec^ 



Vic-l)/c 
2 

/c 
2 



{x — ct) 
(x-t/{l + c)) 



if 



if 



c > 1 

< c < 1. 



(2.6) 
(2.7) 

(2.8) 
(2.9) 



In the (x, t) space the first set of solutions will be solitary wave like and will move 
with a translational velocity c > 1; we call these fast waves. The second set of 
solutions will be periodic, and will have a translational velocity c G (0, 1); we will 

call these the slow waves that are unbounded if ^/c(^x — t/{l + c)j = (2j + 1)tt, 

jez. 

At the boundary between the fast and slow waves we will encounter periodic 
solutions that will be given in terms of elliptic functions that will travel with 
velocity c = 1 as we will see next. 



When c = 1, (2.3) becomes 



+ 0CC = 0. 



(2.10) 



which is not by all means a simpler equation. By multiplying by (f)(^ and integrating 
once again we obtain 

f + (<^c)' = ^, (2.11) 
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where A is some nonzero constant of integration. Now let's use the substitution 
/i^ = A — (/), where /i = m(C)- Hence (2.11 ) becomes 



Moreover, let's assume /i(C) = 3i/'*Ai/2z(C), then (pl2l) becomes 



This differential equation will be solved using Jacobian elliptic functions. 



Lemma 2.1. Ij Z =\ + 2z^ cos 2q; + z'^ , then 



, r dz I , 2x 

u(x,k) — / —= — -sn 

^ ' ' Jo VZ 2 



with k = sin a 

Proof. Putting z — tan 9, we find 
u{x, k) : 



^6* 



'0 a/i - sin^ a sin^ 20 

followed hy y = sin 20, which in turn leads to 

u{x, k) ' ^ ^ 



lo v/(i-y')(i-fc'y')' 

where the integrand is the elliptic integral of the first kind [S]. 



Therefore, solving (2.131, we obtain the solution in z 

^-±sn{A'/y3'^%k), 



2z 



where k — sin 



Stt _ 73+1 



12 2V2 

It follows that in fj, we will have 
2a/x 



1+Z2 

is the modulus of the Jacobian elliptic function. 



±sn(3"i/2^C,fc), 



where a = S^^^A^^^. By solving the quadratic, we obtain 

l±cn(3-i/2aC,fc) 



/i = ±a- 



sn(3-V2aC,fc) 



(2.12) 



(2.13) 



□ 



(2.14) 



(2.15) 



(2.16) 



Hence, the analytical solution to the (1.4) without viscosity — and trav- 
eling wave velocity c = 1 is finally 



L V sn(Ai/2/3i/4(j,_t)^fej J 

l±cn(Ai/2/3i/4(2,_t)_fc) 



A 



(2.17) 
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These critical solutions are unbounded along the parallel lines in the (x, t) plane 
where x~t ~ 77^^, n> 0, and 2K periodic with K given by the complete elliptic 
integral 

- = / ^ = 2.76806 (2.18) 

2 Jo ^Jz^ - ^/3z2 + 1 



The expressions (2. 8), (2. 9) and (2.17), describe the whole class of solitary 
wave solutions, with spectrum c G (0,00); but we shall generally use the symbol 
(j) to mean some specific member of the class, referring to as the solitary wave 
or periodic solution, which is in concordance with [?]. They present periodic and 
solitary wave solutions to ( 1.4 1, with no restrictions on the constant A, which may 



be advantageous when adapting results to physical problems. The validity of the 
equation as an approximate model for water waves in real systems depends on 
the magnitude |u| of solutions being every where small, and so it is warranted 
to require that the velocity c should be chosen such that the waves will travel at 
a speed that is asymptotically close to the linear shallow water speed about the 
point a; = i, or equivalently x* = yJgHt* . Hence, c can not be much larger than 
one. 

The traveling waves solutions obtained by (2.8), (fast waves) are plotted using 



MATLAB function surf and shown in Fig[T To plot the periodic solution given 



by (2.9) (2.17) is more involved. We used MATLAB function ellipj to compute 
the elliptic functions cn and sn of (2.17). These solutions have singularities along 

^/c.{x — t/{l + c)^ = (2j + l)7r and x — t — respectively, and therefore we 

plot them in different domains see Figs (2) [3] 




Figure 1. Exact solutions of the (1.4) given by (2.8) with c = 1.025. 



Solutions to the BBM equation with viscosity 



7 




Figure 2. Exact solutions of the (1.4) given by (2.9) with c = 0.1 




Figure 3. Exact solutions of (1.4) given by (2.17) with c = 1 



2.2. Viscous layers present, i/ > 

For the case when viscous boundary layers are present, we will use an algebraic 



method to solve (1.6). We will use at last two transformations on both the inde- 



pendent variable and the dependent variable on the Abel's equation 



r= 1- 



which will allow us to write (2.5) as 



XXr-X = 2(r-2r5 -l-l-a + ar 5) 
According to [E' the last equation has the closed form solution 



a + (1 — a)T3 



2r 



3(r2 — l)(r — T2 — a) 



(^12aT5 + 12(1 - a)T - 8t^ + 3t' 



(2.19) 
(2.20) 

(2.21) 
(2.22) 



Now, the idea is to use (2.19), (2.20), and (2.22) to obtain an algebraic equation 



in one variable only. Two achieve this, let us equate (2.20) and (2.22) and use the 
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same variable 



in all three mentioned equations. Therefore, we obtain 



a+(l-a)(l-a)-2(l-a;)'' + (l-a;)'' _ 
-3x[(l-a;)2-(l-a;)-a] ~ 

-2(a+bip+x(l-x)) 

(l-a:)[12a+12(l-a)(l-a:)-8(l-a;)2+3(l-3:)3] 



(2.23) 



After few algebraic manipulations (2.23) can be written in the quartic form 

3a;'* - 4x^ - Uax^ + 6(2a - l)x + 6a + 7 + Gbc/) = 
Now, we will use the fact that x — ^^^^ in ( 2.24[ ) and we obtain 

a;'*(0) + a3{(t))uj^{(t)) + a2{(j))uj'^ {(f)) + ai(0)a;(0) + ao(0) 
where the coefficients of the quartic are 

03(0) 



0, 



02(0) 
ai(0) 

ao(0) 



ic 

402(1 - c) 



(2.24) 



(2.25) 



(2.26) 



2v(j)^ {v'^ + 2c{l - c)) 

1/204 |^7z/2--3c(0 + 2(1 -c))j 
3^4 



Since 



(2.27) 



we obtain the solution 0i 
for 0(C) with a;(0) as a root of the algebraic equation (2.251 



ct) implicitly by solving the first order ODE (2.27) 



First, we found the roots a;i234((/>) using MATLAB function solve. The ex- 
pressions are highly nonlinear, extremely lengthy and therefore not presented here. 



As such, the solutions of the ODE (2.27) for the nonzero values of v take very long 



time to compute and due to the possible subtractive cancellation the approximate 
solutions may not be satisfactory. However, we can compute the roots for the cases 
with 1/ = and 1/ = 1 easily. We solve the ODE \2.21\ using the MATL AB func- 
tion ode45 using one of the quartic roots uji for some i of equation (2.25 1 for both 
viscosity parameters p = and 1^ = 1, and traveling velocity speed c — 1.025 and 
c — 0.1. The numerical solutions '/>(C) are plotted as (j){x ~ ct) versus t, and they 
qualitatively agree well with the closed form solutions presented in §2.1. 

In the absence of viscosity v = 0, we obtained both traveling wave solution 
for c > 1 in Fig|4j and the periodic solution for c > 1 in Figj5j In the presence of 
viscosity, 1/ = 1 we expect some damping effect as time progresses. This is depicted 
in Figs|6[7l 

From a qualitative point of view Fig|4] is to be compared to Fig{T] and Fig|5] 
to FigH 
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in using niM). parameters: c-1 .025 and v-D 



Figure 4. Numerical solutions to (1.4) (Algebraic Method) pre- 
dicting the traveling waves solutions for c > 1. 



solution using mM, parameters: csO. 




Figure 5. Numerical solutions to (1.4) (Algebraic Method) pre- 



dicting the periodic solutions for c < 1. 




Figure 6. Numerical solutions to (1.6) (Algebraic Method) pre- 
dicting the traveling waves solutions for c > 1. 
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ig lo^W], paramelers: c-0.1 andv^ 



Figure 7. Numerical solutions to (1.61 (Algebraic Method) pre- 
dicting the periodic solutions for c < 1. 



3. Stability of the viscous waves 



In this section we consider the two-mode dynamical system of (2.3) 



(3.1) 



2^ 



with equilibrium points in the phase plane ((/>, ■)/)) at (0, 0) and 
(— ^,0) — (2(c — 1),0). All the equilibrium points lie only in the plane V' = 0, 
in the {(j),ip,c) space. The bifurcation curves are given by 0(0 — 2(c — 1)) = 0, 
which are two straight lines intersecting at c — 1. There is a bifurcation point at 
c — 1, since the number of equilibrium points changes from two (0 < c < 1), to 
one (c = 1) and back to two (c > 1), as c increases. 
Near the origin (^QjV'o) = (0;0) 

(3.2) 

Following standard methods of phase-plane analysis the characteristic polynomial 

(3.3) 



" ■ 




■ 1 ■ 








1 ■ 




■ " 






7 a 
2 2 . 








1 — c z-* 
c c . 







of the Jacobian matrix of (3.1) evaluated at the fixed point {(f>o,ipo) is 

5o(A) = -poA + go = 0, 



where Po = ^ , and go = ^ ■ 

Since > 0, c > 0, then pg > 0, hence the origin is unstable . Also, putting 



4qo 



-4e(l- 



-, we have the following cases: 



(i) < c < 1, gives go > 0. If > 2y^ c(l — c) > 0, the origin is unstable node 
see Fig. [sj if < < 2^c{l — c), the origin is unstable spiral, see Fig. [oj 

(ii) c > 1, gives go < A > 0, hence the origin is a saddle point, see Figs. 

[Mill 
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Figure 8. Phase-plane for c = 0.5, v ^l, (0.0) node, (—1,0) saddle 
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Figure 9. Phase-plane for c = 0.5, v = 0.1, (0.0) spiral, (-1, 0) saddle 

Near the secondary fixed point (0i, -01) = (2(c — 1), 0) 

(/i 1 r 1 

V> J [ -i-2/3(c-l) f 









■ 1 ■ 




■ " 








1 — C 

c c . 







(3.4) 



The characteristic polynomial of the Jacobian matrix of (3.11 evaluated at the 
fixed point ((/ii,?/;!) is 

5i(A) = A^-piA + qi =0, (3.5) 

where po = f , and gi = -go = - 

Since pi = PO) the second fixed point is also unstable , and moreover we have 
the cases: 
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(i) < c < 1, gives gi < => A > 0, hence the fixed point is a saddle point, see 
Figs-I^ 

(ii) c > 1, gives qi > 0. If < < 2-\/c(c — 1), the fixed point is unstable spiral, 
Fig. [lljor an unstable node for v > 2^ c{c — 1), Fig. 10 

All the other remaining cases, i.e, both fixed points collide, are degenerate, since 
if c = 1, g{\) = A2 - v\. 




Figure 10. Phase-plane for c = 1.5, i/ = 1, (0.0) saddle, (1,0) node 




Figure 11. Phase-plane for c = 1.5, v = 0.1, (0.0) saddle, (1, 0) spiral 

Note that when there is no viscosity, the system has Hamiltonian H{(l),ip) = 
\ip'^ + j(j>'^ + ^(p^- Therefore, along any phase path H{(f),ip) = constant. In this 
case, po — pi — 0, and hence 
(i) ((/'oiV'o) is a saddle and {(jjijipi) is a center when qi > 0, and 
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(ii) {(jjoj^pa) is a center and is a saddle when qi < 0. 

Therefore, in this case the unstable spirals from Figs. [9j [Tl] (which correspond to 
the case with small viscosity v = 0.1) become centers when there is no viscosity 
1^ = 0. 

This is an example of a transcritical bifurcation where, at the intersection of 
the two bifurcation curves (/> = and = 2(c — 1), the equilibrium changes from 
one curve to the other at the bifurcation point. As c increases through one, the 
saddle point collides with the unstable node, and then remains there whilst the 
unstable node or spiral moves away from (ipi, (jji). 



4. Summary and conclusions 



In this paper a basic theory for the BBM equation ( 1.4 1, and its extension ( 1.6 ) to 



include the dissipation term was shown. When the viscous terms are not present, 



(1.4) has traveling wave solutions that depend critically on the traveling wave 
velocity. When the velocity is sub unitary, it has solitary wave like solutions. If the 
velocity is super unitary, the solutions become periodic and unbounded. At the 



interface between the two cases, when c = 1, (1.4 1 has periodic solutions in terms 
of the elliptic functions. Also, and ad-hoc theory based on the dissipative term 
was presented, in which we have found a set of solutions in terms of an implicit 
function that was solved numerically in MATLAB. There is a good agreement in 
the qualitative behavior of the solutions obtained numerically using the algebraic 
method and the closed form solutions. Based on dynamical systems theory we have 



proved that the solutions of ( 1.6 1 experience a transcritical bifurcation when c — 1, 
where at the intersection of the two bifurcation curves, stable equilibrium changes 
from one curve to the other at the bifurcation point. As the velocity changes, the 
saddle point collides with the node at the origin, and then remains there, while 
the stable node moves away from the origin. 
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